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In this paper we consider the matrix representation for the spin operator in the continuum limit 
of the periodic Ising model. We fix the interval [-L, L] as the target for the continuum limit in the x 
variable with periodic boundary conditions <j{L) = cr(—L) for the spin operator. The scaling limit we 
are interested in is a limit in which the lattice spacing, St, for the Ising model in the interval [— L, L] 
tends to as the temperature T f T c approaches the critical temperature. The relationship between 
the lattice spacing and the temperature is determined by taking the inverse, S^ 1 , to be the correlation 
length at temperature T. The interval [-L, L] will then have length 2L in units of correlation length 
throughout the limiting process. The results for the matrix elements of the Ising spin operator will 
have consequences for the correlations in the continuum limit for the cylinder [-L, L] x R and also for 
the torus [-L,L] x [-M,M]. 

We will not attempt to control the convergence of the scaling limit in this paper. This is partly 
because the convergence issue is much simplified if one uses formulas for the matrix elements on a 
finite periodic lattice that were conjectured by Bugrij and Lisovyy (see Q], [3] and [5J). A proof of 
these formulas has appeared but is quite complicated and [5J. We expect that the technique we 
use in this paper can be extended to deal with the finite periodic lattice to give an alternative proof of 
the Bugrij-Lisovyy conjecture. The results for the continuum limit that we address in this paper were 
already announced in [TU] and we should mention that the principal technique we employ is a Green 
function construction that we learned from the paper of Lisovyy [5] . 

We begin by recounting some results from Grethe Hystad's thesis [7J, that provide a representation 
for the continuum limit of the Ising model (under the hypothesis that the Bugrij-Lisovyy conjecture is 
correct). The framework for her thesis is a reworking of Bruria Kaufmann's 1948 paper on the periodic 
Ising model [4] We can avoid some extraneous detail if we limit our considerations at the start to a 
finite lattice A^ m = {—£, ■ ■ ■ , £} x {— m, • • • , m} where I and m are positive integers. A configuration 
of spins on the lattice is a map, 

a : A^ m {-1, 1}. 

In this paper we are exclusively interested in the periodic boundary conditions a(—£,j) = a(£,j) for all 
j. We will be interested both in the boundary conditions for periodic behavior in the vertical direction, 
<j{j,—m) = a(j,m) for all j and also in the cylindrical limit m — > oo. The reader might note that 
little appears to be gained in the periodic situation by having the lattice sites run from — £ to £ rather 
than from to £ as is perhaps customary. However, there are some differences in our treatment of 
the model for odd and even numbers of horizontal lattice sites that make it simpler to confine our 
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attention to the odd case. The energy of a periodic configuration for the Ising model is, 

Ei. m (p) = - Jija(i)a(j) 

(i,3) 

where the sum is over nearest neighbors i and j in Af j?n and in our considerations there are only two 
interaction strengths, Jij = J\ > when i and j are horizontal neighbors and Ji,j — J2 > when i and 
j are vertical neighbors. Of course, in the periodic case lattice points with the same second coordinate 
and first coordinates £ and —I are nearest neighbors as are points with the same first coordinate and 
second coordinates to and —to. The partition function at temperature T is the sum of the Boltzmann 
weights associated to a configuration, 

Zl,m = y^^exp(-E e . m (a)/k B T). 

a 

The Boltzmann constant, ks, appears in this formula, for reasons of tradition, but nothing is lost for 
us replacing Jj by Jj/ks and setting ks = 1 • 

Kaufmann's basic result is a formula for this partition function as the trace of the 2m + 1 power of 
a transfer matrix, V, that has a characterization in the spin representation of the orthogonal group, 

Zt, m = Tr(V 2m+1 ). 



1 Transfer Matrix 

Our first goal is to summarize the reformulation of Kaufmann's result for V that can be found in the 
dissertation [7]. The transfer matrix V acts on the tensor product, 

l 

n = (g) cl 

n=-t 

where is just a copy of C 2 . Suppose that X is a map on C 2 . Let X n denote the linear 
transformation on H that acts as X on the n th slot and the identity in the remaining slots. There is 
a representation of the Clifford algebra on H determined by the action of generators, 




(1) 



where, 

Then {qk,Pk} satisfy the usual generator relations for the Clifford algebra, 

PkPi +PiPk = 2S kh q k qi +qiq k = 2Si k ,qkPi + Piqk = 0. 
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We introduce the vector space W of complex linear combinations of q k and pk with coordinates, 

1 1 

W 3 w = —= ^ x k(w)qk + Vk(w)p k , 
v 2 k=-t 

and distinguished complex bilinear form, 

l 

(w,w')= ^ Xk(w)x k (w') +y k (w)y k (w'). 

k=-£ 

The {q k ,Pk} are then generators of an irreducible representation of the Clifford algebra Cliff (W) [11]. 
The conjugation, 

1 1 

w = —j= ^kQk+VkPk, 

determines an Hermitian inner product on W, 

(w, 11)') — (w, w'). 

The vectors qk,Pk are real with respect to this conjugation and since they are self-adjoint with respect 
to the standard inner product on the tensor product of copies of the Hermitian inner product space 
C 2 this representation of Cliff (W) is a * — representation. 

An important role in Kaufmann's analysis is played by the parity operator, 

U = \\X k = ipkqk- 

k k 

Evidently U 2 — 1 and for reasons that will be apparent shortly we write Ha for the +1 eigenspace of 
U and Hp for the —1 eigenspace of U. The transfer matrix V respects the direct sum decomposition 
H = Ha ®Hp and we write, 

V = V A ®V P . 

In order to characterize the maps Va and Vp we introduce the finite Fourier transforms Fa and Tp. 
Let, 

l e = {-£,-£+ 

and note that in the {qk,Pk\ basis we can think of W as the finite sequence space, 
For / € e 2 (l e , C 2 ) and z G C write, 

Let S A = {z £ C\z 2e+1 = -1} and T, P = {z G C\z 2e+1 = I}. We refer to T, A and T, P as the 
anti-periodic and periodic spectrum on the unit circle. Then J- aJ and J-pf are respectively defined 
as the restrictions of J-f to the anti-periodic and periodic spectral points on the unit circle, 

F A f(z)=Ff(z)forze-£ A , 

Jp/(z)=J/(2)forzeE P . (2) 
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It is easy to confirm the inversion formulas, 

f(k) = , 1 V F A f{z)z- k iork 6 1, 

f(k) = -^= J2 Fp f(z)z~ k for keZt 



For \z\ = 1 define a 2 x 2 matrix, 

T z (F) = e -^o + (z) + e ^o_( z ),withg ± (z) = i(^ ^ Tu | (z) ) 

where 7(2) and w(z) are defined by, 

ch.7(z) = c* 2 c\ — s* 2 s\(z + z _1 )/2 
. . ^i(«)^4 2 («) 

Note that for brevity we write, 

chx = coshx, andsha; = sinha;. 
The constants Sj and Cj and their "duals" s* and c* are defined by, 

Sj = sh{2J J /k B T),c j = ch{2Jj/k B T), 

* — 1 * —1 

s j ~ s j '°j ~ c i s 3 ' 

and finally, 

Aj{z) = y/aj - z 

with, 

ol\ = (cj - s*)(c 2 + s 2 ),a 2 = (cj + s*)(c2 + s 2 ) 

It will be simplest for us to confine our attention to what happens when the temperature, T, is strictly 
less than the critical temperature (T < T c ) 8 ]. In this case a 2 > eti > 1 and smooth square roots 
S 1 3 z — > Aj(z) are uniquely determined by the normalization, Aj(l) > 0. 

Let Ta(V) denote the operator on W whose action in the JFa representation is given by, 

T A f{z) -> T z {V)T A f{z), 

with a completely analogous definition for Tp(V). Let denote the operator on W whose action in 
the Ta representation is given by, 

TaS(z) -> Q±{z)T A f{z), 

with an analogous definition for Qp. 

Both Ta[V) and Tp(V) have positive real spectrum and neither has 1 as an eigenvalue (for finite £). 
Let W\ = Q A W denote the span of the eigenvectors for T A (V) that have eigenvalues between and 1. 
and let W A = Q A W denote the span of the eigenvectors for T A (V) that have eigenvalues greater than 
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1. Let T4 denote the restriction of Ta(V) to the subspace W\. Make precisely analogous definitions 
for Wp and T+ . 

In the dissertation [7] it is proved that Ha is unitarily equivalent to the even tensor algebra over 
W\. That is, 

Ua ^ A\t CVCD (WX) = C © W+ A W% © • • • © A MW * , 
where A. k W A is the k fold alternating tensor product of W\ with itself, and that in this representation, 

V A = \ A (l®Ti®T+®---®(T+)® 2e ), (3) 

where Xa is the largest eigenvalue of Va given by, 

zESa 

In a similar fashion Hp is unitarily equivalent to the even tensor algebra over Wp , 

Hp ~ Alt cvcn (I¥+) = C W+ A W+ • • • © A 2£ HA+, (4) 

with, 

V^p = A P (1 © T+ ® T+ © • • • © (T+) g)2f ), 
where Ap is the largest eigenvalue of Vp given by, 

Ap = cxpi ^ 7(2). 



Remark 1 This representation of T Z (V) differs from the representation in [7] by conjugation by, 




representing translation by 1 lattice unit in the qk basis elements (with qt+\ equal to— q-i, or q-£ 
depending on which transform, Ta or J-p, is relevant). 

This is convenient since the corresponding "conjugation" acting on the induced rotation for the 
spin operator, reduces that operator to the "difference" of translations acting in the periodic and 
anti-periodic sectors. 

2 Spin operator 

The spin operator at (j, 0), which we write as aj acts on the Clifford generators, 

<7jqk<?j X = - sgn(/c - j - l)q k (5) 
UjPkaj 1 = - sgn(/c - j - l)p h 
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Again, compared to [7] this is conjugated by ( ^ 1 ) S ^ n operators <Xj anti-commute with U 

and hence map Ha into Hp and vice-versa. Thus we write, 







r AP 



(Jj = y ^ PA 3 q J acting on Ha © Hp. 

One way to make use of {fi} is to note that H = Ha © Hp is also unitarily equivalent to both, 
Alt(W^) and Alt (Wp). Thus we can regard Oj as a map, 

o j : Mt(WX) -> Alt(Wj), (6) 

each space carrying an irreducible *— representation of the Clifford algebra Cliff (W). The unitary 
map, Ua, that allows us to (projectively) identify H = (C 2 )® 2t+1 with Alt(Wj{") intertwines the action 
of the Clifford algebra on H determined by the generators qt,Pk with the Fock representation Fa on 
Alt(Wjj"). That is, for w — x k<lk + UkPk we have, 

U A w = F A (w)U A . (7) 

Recall that the Fock representation is given by, 

Fa{w) = c{w\) + a(w^), where = Q^w, 

and c(-) and a(-) are creation and annihilation operators [8]. Since are self-adjoint with respect to 
the Hermitian inner product on W, Fa determines a *— representation of Cliff (W), and since W is finite 
dimensional all such irreducible representations are unitarily equivalent. Thus XJ a '■ H — > Alt(W^") 
exists and is determined up to a multiple of absolute value 1 by (0. 

The map Oj is then characterized up to a scalar multiple by the relations (j6]) and we can recover 
a^ A by restricting ([6]) to the even subspace, AltevenCW - ^" ). This is important for us since we are 
interested in understanding the matrix elements of <jj in a basis of eigenvectors for the transfer matrix 
V and these are simple only in the spaces, A\t cvcn (W^ p ) . In a similar fashion we can regard, 

a 3 : Alt(W/) -> Alt(W2), 

and understand a AP as the restriction of Oj to the even subspace AltovcnCW^p )■ The reader should 
keep in mind, however, that this map is not to be confused with ([6]) - in fact, since a 2 — 1, this map 
can be identified with aj 1 for <jj coming from ©. 

A fairly dramatic simplification of the induced rotation for aj occurs for j = t. In this case, ([6]) 
implies that the induced rotation for ai is the identity. Thus ai is an intertwining map for the Fock 
representations Fa and Fp. That is, 

a e F A (w) = F P {w)a L (8) 

One can recover the spin operators Uk by conjugating with the appropriate "space" translations acting 
on, 

Alteven(W^) © Alt ovcn (W+). 

Translation by I lattice unit is given by, 

1 © (z (g) z) © • • • © z® 21 

in both even tensor algebras. However, since z G in the first and z € Sp in the second, the 
spectrum of space translations is different in each summand. This is reflected by the induced rotation 
being 21 + lanti-periodic in W = £ 2 (I m: C 2 ) in the first summand and 21+1 periodic in the second 
summand. This difference produces the signs in (|6j). 
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3 Eigenvectors of the Transfer Matrix 

The operators Ta(V) and Tp(V) acting on W arc sclf-adjoint. We introduce an orthonormal basis of 
eigenvectors for these maps. Define, 



a(z) = 



I Mz)Mz) 



which we normalize so that it is positive for z = 1 (and as usual for us, T < T c so the square roots 
are all holomorphic functions in a neighborhood of z € S 1 ). Then for z G XU> has eigenvalues 
exp(— j(z)) and exp(7(z)) with associated eigenvectors e A (z) and e A (z), given as functions in the Ta 
representation of W by, 



u^e+(z,u) = -t= ^ JJ^-i ) 8(z,u),foiz,ue Ea, 
V2 V -* a ( 2 )" 



u -> e A (^,M) = -^= ( ) 5(z,u), for z,u £ Y,a- (9) 



Here, 



, . f 1 if u = z 
0(z, u) = < f . 

v ' y o if u i= z 

Then {e A (z),e A (z)} ze j: A is an orthonormal basis for W with respect the Hermitian inner product and, 

(e+(z),e A (z'))=S(z,z'), 

so {e^(^)} 2 gSA an d { e ^_{ z )}zeT, A are dual basis for W A and with respect to the complex linear 
pairing between these two subspaces. We make exactly the same definitions for ep(z) except that 
zeEp in this case. 

We are interested in the matrix elements for the spin operator a = an that connect the eigenvectors 
for V (all wedge products of the vectors e A (z) and ep~(z)). So, for example, let 

z = (zi , Z2, ■ ■ ■ , Zk) with Zi € XU and z t ^ Zj for i ^ j; 

and write z £ XU in such circumstances. 
If k is even, the vector e A (z) defined by, 

40) = e + A { Zl )A---Ae + A {z k ), 

is an eigenvector for the transfer matrix V with eigenvalue, 



cxp \ \ 7(z)-5Z 7 ( z )j 

V zeT, A zez J 



We write z e z iff z = Zj for some j = 1, 2, . . . , k, and we write k = Exchanging the roles of P 

and A let, 

= (4>4, • • • >4) with ^ e Spandz- ^ z^-fori 7^ j. 

Then if fc is even, 

e+(z') = e+(4)A...Ae+(4) 
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is an eigenvector for the transfer matrix V with eigenvalue, 



exp 



[\ E 7« - E too) 



4 Matrix of the spin operator 

The matrix element, 

(e%(z'),ae\(z)) for z' e Epandz e E^, 
makes sense if we think of ct as a map, 

cr : Alt(W^) -> Alt(Wp"). 

Of course, e~p(z') and e.\{z) are eigenvectors for V only if both and #z are even. However it is 
easier to explain the version of Wick's theorem that allows us to reduce the calculation of these matrix 
to a few basic types if we admit matrix elements for arbitrary wedge products. 

Let Oa and Op denote the unit vacuum vectors in Alt(W"f) and Alt (Wp ) (the reader might note that 
these vectors are only defined up a multiplier of absolute value 1 by the abstract unitary equivalence 
of these spaces with the original tensor product - we will have more to say about this later). The 
induced rotation T(a) is the identity on W. We write, 



T(a) 



A B 
C D 



for the matrix of this map taking one from the W\ © W A splitting of W to the Wp © W P splitting. 
So, for example D = Qp\ w - , maps W A into Wp. We are particularly interested in the case where D 
is invertible. This happens precisely when (OpjCrO^) 7^ 0. In this circumstance one has the following 
formulas for the simplest matrix elements of a [7] , 



and 



(0p,a0 A ) 

Here, D T : Wp — > W\ is the transpose of D with respect to the bilinear form (•, •) and 



D- T : W+ -> W+, 



is the inverse (D T )~ 1 . 

The following result is an extension of a well known reduction formula for the matrix elements of 
an element in the spin representation of the orthogonal group [8 . 



Theorem 1 Let z E Ha and z' G £p. Then, supposing {0p,(j0a} ^ 0, 



(OpjCTOa) V Rzxz' Rzxz J 



where Pf is £/ie Pfaffian and the skew symmetric matrix R has matrix elements, 



{Rz'xz') 



(e+(z;)Ae+(z;-),<7(U) 
(0 P ,a0 A ) 



-l 



{Rz'xz) 




(ep~(z<),ae A ~{z j )) 
(0p,a0 A ) 




and 




(0p,ae+(^) Ae+(zj) 
(0p,ct0a) 



This result is true even if #z' and #z are not even, but the matrix elements for a no longer connect 
eigenvectors for the transfer matrix in that case. It is evident that to find explicit formulas for the 
spin matrix elements by this method one needs to calculate D . We turn to the principal result in 
this paper, a calculation of D^ 1 in the periodic scaling limit. 

5 The periodic scaling limit 

As mentioned above the limit we are interested in is a continuum limit for a lattice theory that "lives" 
on the real domain [— L, L] x R, with a lattice spacing that tends to as the temperature tends 
to the critical temperature. Periodic boundary conditions on [— L, L] are maintained for the spins 
throughout. For the cylindrical case {M — > oo) the correlations can be understood as 0a vacuum 
expectations of products of spin operators [7] . The doubly periodic case the spin variable lives on the 
torus [— L, L] x [— M, M)/ ~(the equivalence identifies L ~ —L and M ~ — M in the first and second 
factors respectively) and correlations are expressed as traces of the same products. In both cases 
one obtains explicit formulas that depend only on the formulas for the matrix elements of the spin 
operators. We examine what happens to the transfer matrix and the spin operator in this continuum 
limit without worrying too much about convergence questions. The thesis 7J contains results for the 
convergence problem assuming the Bugrij-Lisovyy conjecture [5]. 

The horizontal and vertical correlation lengths are known for the infinite volume Ising model. They 
are the reciprocals of the horizontal and vertical masses defined by, 



The limit we are interested in has Jj/kp fixed and T f T c , where the critical temperature, T c , is defined 



mi(T) = 2K% — 2K*, 
m 2 {T) = 1Kx-2Kl, 



where Kj = Jj/ksT, and K* is defined by, 



sh{2Kj)sh{2K*) = 1. 



by. 
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The critical temperature is defined by the self-dual condition, K\ = K 2 , and one sees that both m\{T) 
and rri2{T) tend to zero as T f T c . We look at the model with horizontal lattice spacing mi(T) and 
vertical lattice spacing m 2 (T). More specifically fix L > and let t(T) — [Lmi(T) -1 ] where [x] is the 
greatest integer less than x. Consider the finite Fourier transform on, 

W = l 2 {I l{T) ,C 2 ), 

written out in terms of the scaled variable x <E mi(T)lM<r) C [-L, L], 

Fffk) = 1 Y f(x/m x (T)) exp ( %— \ . 

The map, 

f ^ F(x) = -^=f(x/ mi (T)), 
V m n T ) 

is a unitary map from the I 2 space on Imt) to the £ 2 sequence space on mi (T)2Wy) , with points 
weighted by the mass m\(T). The naive limit of this Fourier transform as T f T c is, 

FF{k) = -^L [ L F(x) exp f l Ip\ dXi ( 10 ) 



2L 

where k € Z in the case of the periodic Fourier transform and k G Z + 1/2 for the anti-periodic Fourier 
transform. Note that we have changed the look of the finite Fourier transform © by the substitution, 

(2irik \ 
~ 2£ + 1 J w ^ crc ^ e (periodic) or fc e Z + l/2(anti — periodic). 

The parametrization by —A: rather than k is motivated by the desire to get the usual shape for the 
Fourier series coefficients (fTU|) . The scaled vertical coordinate is y € rri2(T)Z. To understand the 
scaling limit we examine the "infinitesimal" generator jQ of T(V). We are interested in the limit, 



lim m 2 (T) 1 j (e m^m^ Q (e wm+i j 



TfT, 

where /eg Z or k E .Z + 1/2 for the periodic or anti-periodic sectors respectively. In [5] it is shown 
that, 

lim m 2 (r) _1 7 ( e" 2 ^ T )+ 1 ) = y/l +p 2 , where p = — . 

Recall that, 

w{z) 
w(z) 

A 1 (z)A 2 (z) 



Q(z) = - 

with, 

w(z) = i- 

A 1 {z- 1 )A 2 {z- 1 ) 

Since a 2 > 1 in the limit T | r c we see that, 



-4 2 (e 



g 2/(T) + l 

lim 



TtTc A 2 (e 2 *w+A 
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However, a% — e rai ' T ' tends to 1 as T f T c and we have the asymptotics, 



Ai «"i ^ ^ ^/mT(r)\/l + where p = 

Thus, 

1 V / \/l+ip . irk 

lim t — — -. — r— = , , tor p = — . 

TtT c Ai L^mTfT ) y/l-ip L 



The induced rotation for the transfer matrix thus scales to the relation between infinitesimal vertical 
translation and infinitesimal horizontal translation, 



d ( 1 \ / -i \ ( -i \ d / 



dx 2 V 1 J p ' V 1 J V -* 7 fei ' V 1 



where we identify ip with d/dxi based on (|10p . Consider the differential equation, 

^ = 0. (11) 



# , fo iU / o i 

dx 2 I i U ) dx ^ \ -i 



After making the substitution, 



this differential equation becomes, 



i _ -2a 

-25 1 



V = 0, (12) 



where, 

2 Vctei 9x2/ ' 2 \dx\ 8x2 

Equation f|12|) is a convenient form for the (Euclidean) Dirac equation in two dimensions. In fact, it 
is useful to generalize this just slightly. We introduce a mass m > into the Dirac equation, 



m —2d 
-2d m 



tp = (13) 



This arises naturally if one introduces the scaling variable mx\ instead of X\. The spin correlations 
depend only on the product mL so L and m are not really independent parameters. However, in 
the planar case the short distance behavior of the scaled correlations can be studied by analyzing the 
m — > behavior of solutions to (|13p and we expect something similar is possible in this case. 
In the scaling limit it is natural to parametrize periodic and anti-periodic spectrum, 



L 

and 



{^n:nez}, (14) 
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Of course, this notation conflicts with the earlier definition of S p t A but for the remainder of the paper 
we work with the scaling limit so the possibility of confusion is small. The scaling limit of the space 
W\ can be identified as the space of solutions, ip(x) — ip{x\ 1 X2), to (fl~3|) that are 2L anti-periodic in 
x\, are defined in the upper half plane ^s{x) = xi > with L 2 [— L, L] boundary values for X2 = and 
tend to as X2 — > +oo. is similar but the solutions are defined for X2 < and tend to zero as 
X2 —> — OO. The limits for W$ are similar but are defined in terms of 2L periodic solutions for (fT5|) . 
We understand aj, as an interwining map, 

a L : Mt{WX) -> Alt(Wj). 

This needs some further elaboration since in the scaling limit the Fock representations, Fa and Fp 
are not unitarily equivalent. Nonetheless, for the purposes of finding matrix elements of the molified 
spin operator V e aj_y e (for e, e' > 0) and formulas for the spin correlations it suffices to invert the D 
matrix element of the identity map from W\ © to Wp © Wp . This is the problem to which we 
now turn. 



1 
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6 The spectral curve for the Dirac equation 

A crucial ingredient in our analysis of the scaling limit is a spectral curve for the Dirac equation that 
allows one to analytically continue the eigenfunctions that span W\ P to the eigenfunctions that span 

For u £ C the function 

e(x, u) := e -f 
is an exponential solution to the Dirac equation 

Ve= { -25 m J e = ' 

where x = x\ + ix% and d = d x = h{d/dx\ — idjdx^)- It is convenient to make the substitution 
u = —er s and to regard the complex strip 3s € [— 7r, n] with the upper and lower edges identified, 
s + 2ni ~ s, as the "spectral curve" for the Dirac equation, Ve = 0. We write, 

£(£>) = C/2-kiZ, 

for this periodic strip. 

For this parametrization the exponential solutions, 



E(x, s) — e(x, —e s ) = e 



e 



that are purely oscillatory in the X\ variable arise for 3s = This fact is, in part, the reason we 

make the choice u — — e~ s . Define 

I 7 7T 

M± = <sGC:^s = ± — 
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Fourier analysis allows us to synthesize the solutions to the Dirac equation of interest to us as linear 
combinations of exponential solutions, E(x 7 s), for s <G M±. Write 



ITT 

s± = s± y 



Then it is easy to check that, 



E{x, s±) — exp(mch(s±)a;i + imsh(s±)x 2 ) y e ~ s ± J = cxp(±imsh(s)xi =p mch(s)x2) y z^i e - 

For s real we see that E(x, s + ) is exponentially small at infinity for x 2 > and E(x, s_) is exponentially 
small at infinity for X2 < 0. It is useful to record a related result for the norm of E(x,s). Let 
x = x\ + ix2 and s = u + iv denote the splitting of x, and s into real and imaginary parts. Then after 
a short calculation one finds, 



||i?(x, s)|| = yl + e 2u exp(ma;i ch(u) cos{v) — mx 2 ch(u) sm(v)). (16) 



7 Periodic and anti-periodic Green functions 

Next we calculate the Green function for the operator D acting on the smooth C 2 valued functions 
on, 

[-L, L] x R, 

with either periodic or antiperiodic boundary conditions on [— L, L\. Introduce the Fourier coefficients, 

/(P,6 = ^ J dx! J dx 2 f(x 1 ,x 2 )e- i ^e- i ^\ 
The inversion formula for smooth periodic functions in x\ which are small at oo in X2 is, 



/CO 
-oo 



The inversion formula for smooth anti-periodic functions in x\ which are small at oo in X2 is, 



/oo 



The Green function for T> with a domain that contains the smooth periodic functions of x\ G [— L, L] 
we denote by Gp. The Green function for V with a domain that contains the smooth anti-periodic 
functions of x\ e [— L, L] we denote by G a- For x = x\ + ix 2 one finds, 



m 2 + p 2 + £ 2 \ *P _ £ TO 



If x 2 > the £ integration can be "closed" in the upper half plane and if x 2 < the £ integration can 
be "closed" in the lower half plane. One finds, 

GpA*) = 1t £ 4t(v "\„ liP + LU{p)) )e^^iorx 2 >0 1 (17) 
4L r V w(p) V HP-w(p)) m J 
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and 

GpA*) = ±t E ~h( it '^^le^lbr^l (18) 

To simplify some calculations with these Green functions we introduce, 



9p,a( x ) = Gp,a(j:) 



i 
-i 



The functions gp.^x — x') have the advantage that their columns are solutions to the homogeneous 
Dirac equation as a function of x and their rows are solutions to the homogeneous Dirac equation as 
functions of x' (for Gp a one would need to introduce homogeneous solutions to the transpose of the 
Dirac equation to describe the rows). 

Make the substitution p — m sh s in equation (|T7f and multiply on the right by ( ^\ q ) • O ne 



finds, 



9pa{x) = ^- Y 4~ ( &S - -s ) e msh( - s)xi - mch( - s)x Horx 2 > 0, (19) 
AL *—f chs \ —i e s J 



where we've written, 

s G Sp A for roshs G Sp,^. 
Make the substitution p = — mshs in equation (|18JI and multiply on the right by ( ^ . * I. One 



finds 



1 NT- 1 



gpA (x) = — -r-i _ ■ _ e - s ) e -™ lsh ( s )*i+ mch ( s )^for:E 2 < (21) 



Using (Unj), (|2"Tj) and the definition of E(x, s±) above, one can verify, 

g P , A {x-x') = -1 V -f-£(a;,« + )£;(s / ) fl_) T for9f(a;-a/)>0, (22) 

4_L ch s 

and 

g P:A (x-x') = -j I - E ^^O^s+ffor^-a/) < 0, (23) 

where JT T is the transpose of X . These formulas will be useful when we turn to our principal concern- 
a formula for the Green function of the Dirac operator on the strip [-L, L] x JSwith periodic boundary 
conditions in the lower half plane (x2 < 0) and anti-periodic boundary conditions in the upper half 
plane. 

We introduce periodic and anti-periodic spectral transforms that are closely related to the repre- 
sentations (f2"2"| and (j2"3"|) . For f(x) a smooth periodic (P) or smooth anti-periodic (A) function of 
X G [-L, L], define spectral transforms Sp and <S^ 

r-L 



Si, P f(s) = J L e s / 2 E(x,s.) T f(x)dx fors e V A<P , 
5J,p/(«) = ^ y e s/2 £(z, s+) T /(^)dx for a G £' AP . (24) 
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The Fourier inversion formula implies that, 

= E (S+f(s)E(x,s + )+Spf(s)E(x,s.)) 

f& = E (^/OO^ «+) + «-)) (25) 

Where /(a;) is 2L periodic in the first case, and 2L anti-periodic in the second case (the formulas 
are true in an appropriate L 2 ,£ 2 sense even if the function f(x) is neither periodic or anti-periodic 
but the convergence of the inversion sum is especially good if say the function / is smooth and (anti- 
periodic and the spectral transform is the corresponding (anti-)periodic transform). The formulas 
(123)) respectively determine the splitting of / in the Wp~ © Wp and W\ © WT representations of 
W = L 2 [—L, L]. The scaling limit analogues of the projections Q\ p are, 

Again, since in this part of the paper we are exclusively concerned with the continuum limit there 
should be little danger of a confusion with the earlier usage for p . 
The Plancherel theorem implies, 



'p. 



(26) 



k J_ L i/(*)i 2 ^ - E d(s) (WW 3 + |5 - /(s)|2 ) 

It is useful to introduce the components of the inverse spectral transform, 1Z, as maps, 

K P f(x)= £ J^Lf( s )E(x,8±), (27) 

where / € 1 2 {Y^' A ) or / € £ 2 (Ep) respectively The Plancherel theorem implies that the maps 7\L^ P 
are all isometries from the appropriate square summable sequence space (with weight (2chs) _1 ) into 
L 2 [-L,L}. 

In preparation for the construction of a Green function for the Dirac operator that mixes peri- 
odic and anti-periodic boundary conditions we follow Lisovyy [9] and turn next to the solution of a 
factorization problem on the spectral curve. 

8 Factorization of ctp/cta on the spectral curve 

The function, 

a P {s) = ( e mLchs - e- mLchs )/2 
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has simple O's at the points s — t ± in/ 2 € A4± with 

t e Tf P . 

The function, 

<7 A (s) = ( e mLchs + e - mLchs )/2 
has simple O's at the points s = t ± i7r/2 g M.± with 

Their reciprocals have simple poles that can be used to express sums over the points in the periodic and 
anti-periodic spectrum as contour integrals. Lisovyy's construction of the Green function of interest to 
us depends on the solution of a factorization problem that we describe next. Consider the functions, 



and, 



fo(s) — defined for — e < 3(s) < e, 



/ir(s) = — — — 7-r defined for — ir ^ 3(s) < — n + e or ir — e < 3(s) ^ it 



The precise value of e > is not important but we do want to choose e < it/2 to stay away from the 
zeros of <ja{s) and ap(s). The functions <ja{s) and ap(s) are 2ni periodic in s so we may regard /o 
as a function on £(£>) defined in a neighborhood of S*(x) = and f„ as a function on S(2?) defined in 
a neighborhood of SJ(a;) = ±7r. 

The change of sign in the definition of f v is important because of the particular technique we use to 
produce a holomorphic factorization of the pair (/o, f„) on the spectral curve. As defined /o(s) tends 
to 1 as 5R(s) — > ±oo on the line 3(s) = and /,r(s) tends to 1 as 5ft(s) — > ±oo on the line 9(s) = ±7r. 
In fact, it is useful to note that, 

fo(s) = l + 0(e- michs )for3(s) =0ands->±oo 
/„(« + »tt) = 1 + 0(e- michs )for3(s) = Oands -> ±oo (28) 

As a consequence both log(/o) and log(/ ff ) will tend to as 3?(s) — ► ±oo and this will insure the 
convergence of the integrals we use to define the additive splitting of the logarithms. In fact, since 
neither fo or f v has either a zero or a pole in its simply connected domain of definition they both 
possess holomorphic logarithms which we can normalize to be at oo in s. We henceforth denote 
these choices by log/o and log/^. It follows from (|28p that, 

log/ (s) = 0(e~ michs )for3(s) = Oands ^ ±oo 
log f v (s + in) = 0(e- mLchs )ioT^{s) = Oands ^ ±oo (29) 

We can use the Green function for the Cauchy-Riemann operator on the spectral curve, S(2?), to 
additively split this pair in the usual fashion. Define the level sets for the imaginary part of s, 

Y a = {s|3(s) - a} 

positively oriented by the one form ds. For < S*(x) < ir define, 

F + {z) = ± I ' -fL-iog/,,^- ±- J y -jL_Xo g Us) ds . (30) 



1G 



For -7T < < define, 

F-(z) = ^-[ -^_log/ 7r ( S )d S --L / ^l^bg/ ( s )d S . (31) 

27TZ 7y_ e — e 27TI Jy- e s — e z 

Because all the integrals involved converge absolutely, the functions F±(z) are holomorphic in their 
domains of definition. In fact, since 

1 e s 



2ni e s — e z 

is the kernal of the Green function for the Cauchy-Riemann operator on the spectral curve one can 
see that F + (z) and F_(z) have boundary values such that, 

F+{z) + F_(z) =log/o(z)for3(z) = 

and 

F+(z) + - logf w (z) for 9(z) = ±tt, 

where Y^ is identified with Y-^ in the second equality of boundary values. 

Once one knows this, it is clear that both F+(z) and F_(z) extend holomorphically to the strips 
— e < $s(z) < tt + e and — tt — e < 3(z) < e respectively since the domain of analyticity for log/o and 
log naturally enlarges the domain of analyticity for F± . It is convenient to regard the enlarged 
domains for F± as subsets of the periodic strip E(2?) = C/2iriZ. Thus both F + (z) and F_(z) are 
holomorphic in a neighborhood of both Im(z) = 0, and 3(z) = ±7r in 2(2?). 

Rewrite the formulas (|50|) and (j3"Tj) using s to parametrize Yo and s ± Z7r to parametrize 1± T . One 
finds, 

F ^ = 2r7i ah^ l0 S^> fo r° < < ^ ^ 

and 

F-{z) = ~ !°° 1 bg^da,for-7r<3(s)<0. (33) 
2tt« 7-oo sh(z - s) ovl(s) 

Since the hyperbolic sine is i7r anti-periodic it follows from this representation that, 

F + {z + m) - F_ (z) = 0, for -tt < 3(z) < 0. (34) 

Next we introduce the modifications needed to holomorphically factor <jp jo a with no sign change 
on y±7r. Define, 

X P (s) = e s/2 exp(F + (s)),for-e < 3(a) < tt + e. 

A A (a) = e s/2 exp(-F_(s)),for-7r-£ < 3(a) < s (35) 

We regard both Xp and A a as holomorphic functions defined on open subsets of 2(2?). Keep in mind, 
however, that because e s ^ 2 is 2wi anti-periodic the representation for \p(s) becomes, 

A P (s) = -e s/2 exp(F + (s)),for-7r < 3(a) < -tt + s, 

with a similar modification in the formula for Xa{s) for tt — e < 3(a) ^ tt. 
With this 2tti anti-periodic adjustment one finds that, 

= for s in a neighborhood of Y and Y v = YL n .(inS(D)), 

Aa(s) cr^(s) 
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or more conveniently for us, 

Xp(s) _ Xa(s) 



for s near Y and Y„ = F_ 7r (inS(P)). (36) 



crp(s) cja(s) 

9 The Green function Gp/a{ x , x ') 

Following Lisovyy [S] we write down a Green function for the Dirac operator on [— L, L]x R which has 
periodic boundary conditions on [— L, L] x R + and anti-periodic boundary conditions on [— L, L] x R , 
where R + is the set of positive real numbers and R- is the set of negative real numbers. Define, 

= Y/,— a ~ ^b+a, 

which is the "counterclockwise" oriented boundary of the strip b — a ^ S(z) ^ b + a, centered at Yj, and 
having width 2a. Next we introduce partial definitions for the Green function Gp/a( x : x> ) 01 principal 
interest for us. For ^s(x) > > ^s(x') define, 

Gp/a(x,x')=9a(x-x')-^ I dsj dtE(x,s) Xpis) eS - et t XA fh (x',tf. (37) 

S2ir z J Y a j Y b <7p(s) e s + e l <TA(t) 



For Q(x) < < 3(x') define, 

g p/A ( x , x ')= gp ( x -x')- ds dtE(x,s) AK1 P ) > E(x\t) T . 

32n z J Y * , JY b , aA{s) e s + e f ap(t) 



(38) 



In these formulas we suppose that, 



< b < a < 7r/2. 



In both integrals above the difference ?s(s~t) is ir ± (a — b) or 7r ± (a + b) or — ir ± (a — b) or — 7r ± (a + b) . 
In no case is e s + e* = since this would require that s and £ have imaginary parts that differ by some 
odd multiple of 7r. Also, the restrictions on a and b keep s and t away from the zeros of up and cta- 
Thus the integrals do not have any local singularities. Next we check convergence for large values of 
s and t. It will suffice to illustrate the estimate for the integrand in the s integration in (157)) . Write 
s = u+ iv, and x — x\ + 1x2- Using the easily confirmed observation that F±(s) tends to as \s\ —¥ 00 
we find the asymptotics (u —¥ ±00), 



-E(x, s)\ = O [ -. exp(ma;i ch(u) cos(zj) - mx 2 ch(u) sin(w)) 



'ap(s) ' \ e mLch(u) 

As long as X2 ^ 0, sin(w) ^ 0, and \x\ \ < L, it follows that, for any e such that < e < m{L—x\ cos(i>)), 

,Ap(s) 



a P (s) 



E(x,s)\=0 (exp(-e ch(u))) . 



Since sin(7r/2±a) ^ for < a ^ 7r/2, the s integral in (|38l) converges at 00 provided only that X2 ^ 
and < L (this is true even for the limiting case a = tt/2). In a similar fashion one can see that 
for x' = x' x + ix'2 the t integral in Q38p converges at 00 provided x' 2 ^ and (a;^! < L (again even for 
a = n/2). 
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This is useful since we understand Q37p and (|38[) in two ways. First by taking limits a — > n/2 and 
b — > 7r/2, and then by collapsing the integrals to the residues at the zeros of ap and a a- 

First we take the limit a — > n/2 of the right hand side of ([37)) with b fixed. The analyticity of the 
integrand in s and the asymptotics for large s justify the replacement Y®, 2 — > Yq — Y n in this limit. 
We cannot do the same with the limit b — > n/2 since the points s S Y„ and t — s — in € Yq and s S Y 
and t = s — in € Yl^ are zeros of e s + e*. In the limit b — > 7r/2it is possible to skirt round t = s — in 
with a half circle of radius e. The limit e — > then falls into two pieces. The first is a principal value 
integral and the second is a half residue. The principal value contribution to the limiting value of the 
integral on the right hand side of (|37|) is, 



f , , . \p(s) e s — e f XaU) ^, , xT 

■Km— -j / dsdtE(x,s) P A ff a'.f f 



where, 

r(e) = (y -^) x (r_ w -y )\{(«,t) : it-s + H < £ } 

The half residue contribution is, 



T 



£-/ dsE( X ,s) ^f^f^ E( X ',s-in) 
1D7T Jy -y„ <Jp(s)<7A{S - m) 

Consulting ([M]) we see that, 

Xp(s)Xa(s — in) — —ie s , 

so that this residue contribution is, 

~7?-/ dsE(x,s) 6 — S(x',s - i7r) T . 

This integral collapses to a sum over the residues at s E E P + in/2 and s 6 + in/2. One finds 
(recall that s± — s ± in/2), 



or consulting 

-gp(x-x') - -g A (x-x'). 

Thus for 3(x) > > 3(x'), 

g PM (x,.xO=.gp + A(x-xO-hm-^ / d S ^(x, S )^4^-4^|^',if, (39) 

e^o 32-7T 2 7 r(e) erp(s) e s + e* £ta(*) 

where, 

9p+a = (gp+g A )/2- 
A precisely analogous calculation shows that for 3(x) < < 3(x'), 

S P/A x,x' =5p+a x-a: ') - Km — 5 / dsdtE(x,s) t P )J. E(x',t) T , 

-^0 327T J r / (e) cta(s) e s + e £ crp(t) 
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where r'(e) = (Yl^ — Yq) x (Yq — l / " 7r )\{(s,t) : \t — s — iir\ < e}. We transform the last integral using 
(|36|) ; we replace Xp(t)/ap(t) by \A(t)/<JA(t) and \a(s)/<ja(s) by Ap(s)/crp(s) and then use Y n = F_ w 
on £(£>), to replace r'(e) with T(s). One finds that for 3(x) < < S(x')> 

^^^-^(x-xO-lim-^ / dsdtE(x, s) Xp( f\ ^ ~ e \ XA f E(x', t) T . (40) 

Note that in the formulas (j3T))) and (|40p the integrals on the right hand sides now have the same shape. 
This has important consequences for a polarization that we will associate with the kernel Gp/a( x , x'). 

We turn to the other evaluation of Gp/a( x t x ') that is important for us. We can evaluate the 
integrals in (1571) by residues at the poles on s g Y_/ 2 and t g Y_^/ 2 n °t forgeting the poles in (e s + e*)~ 1 
that arise for s = t + wr (this pole gives rise to — <?a(£ — x '))- Since we want to maintain a > 6, the 
integral residue calculation should be done first. For Q(x) > > SJ(x') we find, 

avxfe-o-sb E E *<-)^^^-~.o E(l <, t r 

ses+ tes A 

where, 

E± A = E PiA ± i7r/2. 

Recalling that s± = s ± ?7r/2, we can turn this into a sum over real arguments, which makes it easier 
to see the convergence of the sum. For ^s(x) > > 3(x'), 

= 52* E E ^M^^^"'^.*-)' (41) 

An analogous residue calculation for G(x) < < SJ(x') shows that, 

^ p M (x,x') = — ^ V V £( IlS )^^M) e - L(ch S+ ch i)jB / / t) T 

P/Ay ' ; 8mL 2 ^ ^ V ' ; shs e s + e* sht V ' y 

ses A tes+ 

And again this transforms into a sum on real arguments. For !J(x) < < Ss^x'), 

= E E <«> 

ses^ tes' p 

The reader can now check that the representations (j39|) and (|4Tj) for Gp/a( x i x ') extend continuously 
from the domain 3(x) > > Sy(x') to 3(x) > ^ 3(x') and the representations (|40[) and (|42|) for 
Gp/a{ x i x ') extend continously from the domain Q(x) < < Sy(x') to the domain 3(x) < ^ Q^x'). 
In each case, we further require that \x%\ < L and \x[\ < L to guarentee that the simple estimates 
we gave above assure pointwise convergence in the integral representations (the sums don't have a 
problem). 

Associated with the Green function Gp/a there is a splitting of L 2 [-L, L] which is important for us. 
We will first discuss this in a naive setting and then introduce the estimates needed to make it work 
mathematically. We hope this will make it easier for the reader to see the reason for the definitions 
dSTJ and (J2SJ). For x, x' e [-L, L] define, 

VpiaK*) = lim / Gp/a( x + fe, x')f(x')dx' (43) 
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and ^ 

V p/ A f( x ) = - lim / Sp/a (x-ie, x')f{x')dx' (44) 

Note that (gTJ) and suggest that, 

€ fl£ , and V p/A f G . (45) 
On the other hand ([39]) and (40]) suggest that for 3?(x) = S(x') = 0, 

Gp/a( x + *0, x') - Gp/ A {x — £0, a:') = ,gp + ,4(a; — x' + iO) — gp+A{x — x' — iO) — 5{x — x 1 ), (46) 
since, 

gp{x — x' + iO) — gp{x — x — iO) — gp{x — x' + iO) — gp(x — x — iO) = 5(x — x'). 
Thus we expect that, 

H/ A f + T > p/Af = f> ( 47 ) 

We will prove later that, 

v% A n, A = o- 

Combined with (|47|) this implies that Vp/ A is a projection. These projections will allow us to explicitly 
invert the map, 

Q P :W A ^W P . (48) 

Note that if, 

' A B 



C D 

is the matrix of the identity map on W from the splitting W A © W A to Wp © Wp then d48]) is just 
the map £>. 

To see how to invert D suppose that / G W A and Qpf = g G W^T. Then f = gp + g where, 
g P G Wp, and so it follows from (1431) . (|47l) and the fact that "Pp^ is a projection that is on Wp 
that, 

/ = T>P/ A f = V p/a G?p + 9) = V P/A g = D~ l g (49) 

Thus (j4T)]) gives an explicit formula for the inversion of (|48|l . 

There are a number of matters that we glossed over in these "pointwise" calculations. The first 
matter we take up is a proof that Vp, A is bounded in L 2 . It is enough to illustrate this for Vp, A using 
the formula (|41jl . Suppose that / is a smooth compactly supported function on [—L,L]. Define, 

Because we are interested in L 2 bounds we can and do ignore the unitary factor e imL ( sh *- shs ) in gj). 
Using the formula ([2"4"| for the spectral tranform and the formulas (|35p for Xp and A a we find, 
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Introduce, 

aw 2mL ^ e s -e* chi pJW v ' 

tez' A 

Then the formula (|2"T|) implies that, 

Because 7?.p is an isometry it will suffice to obtain a suitable estimate for, 



^ 2chs ' 

However, since fs^ft, and e~ F+( ' s +<' _i;, ~( t ~- ) are uniformly bounded in sandi for s G S' P and f G S^, 
and (chi) -1 is summable for t G S' P we can use the Cauchy- Schwartz inequality in ([50)1 to see that, 



Mfl^< C V — V lS r m ' 
^ 2chs ^ ^ chs ^ 2chi 



This finishes the proof that Vp/ A is bounded in L? on the smooth functions of compact support in 
[— L, L], It then extends by density and continuity to a continuous map on all of L 2 . 

Next we show that Vp~/ A f G Wp . Truncating the s sum to a finite range in the contribution that 

the double sum in (|4Tj) makes to the calculation of (j43|) we see that the result is clearly in Wp . Since 
the s sum converges in L 2 we see that ~Ppi A f G Wp. The same argument shows that Vp^ A f G W A . 

One other unresolved issue we wish to consider is the "pointwise" cancellation that produced |46|) . 
We do not have estimates for the terms that cancelled at the endpoints x = ±L and x' = ±L. 
However, the L 2 continuity of the maps Vp^ A makes it natural to define them as L 2 limits defined first 
for functions, /, of compact support on [~L,L]. If / is such a function then the contribution made 
by the terms g A +p in © and gDJ) to Pp~ /A f + V p/A f is, 

\{Qt + Qt) f + \{Qp + Q- a ) f = f- 

Applied to a function of compact support in [-L, L] the issue of convergence for the x' integration 
for the residual kernels in and P0|) in the calculation of Vp/ A j 'does not encounter any difficulties 
at x' = ±L. The resulting functions clearly cancel except possibly at x — ±L. However, since the 
kernels define bounded operators on L 2 the cancelation is good in L 2 and it follows that, 

Finally we take up the proof that Vp/ A maps Wp to 0. Suppose that, 



f(x) = n+F( 



•J' 



where F(s) is non zero for only finitely many s G T/ P . Such / are dense in Wp so it will suffice to 
prove that Vp^ A f — 0. Use (H2l in the definition (PHI) . One encounters the integral, 



L 

E{x',t + ) T f{x')dx'. 

L 
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This can be rewritten as the integral of a one form, 

J Ei(x',t + )fi(x')dx' + E 2 {x ,t + )f 2 {x )dx . 

It is not hard to check that as a consequence of the fact that both x' — > E(x' , t+) and x' — > f(x') are 
well behaved solutions to the Dirac equation in the upper half plane, the integrand in this last integral 
is a closed one form. This means that the contour in the integral can be shifted from $s(x') = to 
^s(x') = R > for any R. However, both solutions tend to as R — >• +oo and it follows that the 
integral itself must vanish. This shows that "Pp^ A f = and finishes the proof. The only change 

needed in the proof that Vp^ A maps W A to is the observation that the product of two anti-periodic 
functions is periodic. This allows the contour deformation argument to proceed. 

Next we turn to an an extension of the definition of Qp/a{ x i x ') to the remaining possible values 
for the argument (x,x'). First we consider the extension of Qp/a( x i x ') to ?${x) > ^ 0. Start 

with (|41[) . In that equation replace the t sum by a contour integral on the boundary of tubular 
neighborhood of 3?(t) = — tt/2 using the kernel a\p(i) _1 . Note that the contour integral that it 
is natural to introduce contains an "extra" pole at the zero of e s + e*. We compensate by adding 
gp(x — x) to the contour integral to obtain a precise representation of the original sum. Suppose that 
^s(x) > and 3(:r/) = 0. Extend the contours in the t integral to the horizontal lines = 0, — tt. 
Replace the ratio A^^/WW in the integrand by \p(t)/ap(t) using equation (plrJ)) . Then collapse 
the resulting integral to the sum of the residues on = tt/2. Keeping in mind that the orientation 
of the contour flips going from = —tt/2 to 3(f) = tt/2 one finds, 

&,* ! 'i=«*-'0-ssE E Ei*,.)*g£**&«-**~*>Bv,tr. (5i) 

se£+ tes+ 

In this form the second term on the right hand side has a continuous extension to ^s(x') ^ 0. A little 
thought shows that the right hand side of (l5Tj) represents the Green function of interest when both x 
and x' are in the upper half plane. An almost identical argument produces the following formula for 
Gp/a( x , x ') when both x and x 1 are in the lower half plane, 

We finish this section by using the projections Vpi A to give formulas for BD^ 1 , and D~ 1 C. Observe 
first that since, 

(c d)> 

is the matrix of the identity, every element / € W A can be written, 

f = Df + Bf, 

or since D is invertible, 

Wa 3 / = 9 + BD- X g, where g G Wp saidBD^g £ Wp~. 
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This is described by saying that W A is the graph of BD 1 : W P — > Wp over W P . Now apply Pp/ A 

to both sides of this last equation and use the fact that W A is in the kernel of this map and that the 
elements in Wp are fixed by this map. One finds, 

= V+ /A g + BD^g for g G W P . 

Thus, 

BD- X g = -V+ /A ghvg € Wp (53) 

Note that (|51|) is a particularly effective representation for Qp/a in the evaluation of this representation 
for BD^ 1 . In particular the gp term in (|51j) makes no contribution since Qp<? = for g € VFp . 

Next we consider a representation for D~ 1 C. Observe first that since the identity map is an 
orthogonal map we have, 

/ A d \ / nr dt \ / i n \ 

(54) 







/ D T 


B T > 












[ C T 




1 = 1 


v 





so that, 



D T B 7 



is the matrix of the identity from W = Wp @ W P to W = W A © W A . Thus for / e Wp we have, 

3 f = D T f + C T f, 

or introducing g = D T f G W A we have, 

^ 9 / = g+ C r L>" T .gfor.g G 
Using the identity CZ? T + Z?C r = that follows from (|54|) this becomes the graph representation, 

W$3f = g- D- l Cgioig G W A . 

Applying Vp^ A to both sides and noting that Wp is in the null space and W A is in the range of this 
projection we find, 

D- 1 Cg = V p/A gioTge W+ (55) 



10 Matrix for D \BD , and/} C in the spectral variables 

In this section we translate the representations (|49l) , ([53)) and (|55|) into the spectral variables ([24|) and 
(125]) . Then we make some adjustments in these representations that we held off because incorporating 
them earlier on would have forced us to keep track of branches for the square root of sh i, or irritating 
factors of exp(i7r/4) in the calculations of the preceeding sections. 

We start with the calculations for D^ 1 . Substitute (|42j) in (l49l) and make obvious use of (|24l) and 
(f2"5]) . One finds the matrix representation of D^ 1 is, 

~ e -m L ^- S/2 Y s) ^ ~ & \ e ~ S 'y {t) e-^Ho.s G E A ,andi G £+ (56) 
2mL e s + e* sh< A p v ' 
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This should be understood in the following sense. Replace s by s_ with the new value of s chosen in 
T,' A . Replace t with t + with the new value of t chosen in Yl p . The matrix of D _1 in (|56[) is implicitly 
indexed by £^ x X P . 

Substituting ([51]) in ([53| one finds (using the convention just explained) that the matrix of BD -1 

BD-- ~ ^ ie -m^ s e-^Xp(s) * - e* e-^A F (t) ^ ^ + 

2mi e s + e* shi p p 

Substitute (|52p in (|55p and one finds for the matrix of D~ 1 C, 



2toL e 5 + e* sh i 



Next we get rid of the asymmetry in these results that comes from using the £ 2 space on the spectral 
side with weight (2cht) . One can pass to the matrix representation in the unweighted I 2 space 
by multiplying on the left by (2chs)~3 and on the right by (2cht)2 where (s,t) £ P x H' A p in 
whichever combination is appropriate. Do this and then adjust the spectral transform further by 
multiplying by e l7r / 4 in the W7 P subspaces and by e~ l7r / 4 in the W Ap subspaccs. One finds, 

Theorem 2 Suppose that 

A B 



C D 

is the matrix of the identity map on W from the W A © W A splitting to the W P © W P splitting. Then 
relative to the modified spectral-inverse spectral transforms, 



^,p/( s ) = ^I-l zmr s E ( x > s -) T f( x ) dx fOTS e z> A .p, 
s A , P f(s) = ^f- L zSr s E ( x > s +) T f( x ) dx fors e E kp- 

and 

f(x) = E seS ' p Z& (e™/ 4 S+f(s)E(x,s + ) + e-™ ^S p f(s)E(x, s_)) 
/(*) = i^t (e^ /4 S+f(s)E(x, s + ) + e^S A f(s)E(x, *_)) 

The matrix elements of D~ x , BD^ 1 and D~ X C are given by, 

~ J- e -m^s. e-'-l^s.) e'- - e'+ e^MU) e - mL *t + fa( f) £ ^ x 
2mL Vch7 e s -+e'+ Vcht V ' A P 

BD~ X ~ - J- e -^+ e- + /^p( a+ ) e- + - e'+ e^)^ for(M) e Ep x jv 
2mL Vch7 e s ++e*+ v/chi p p ' 

2mL Vchs e s -+e*- Vcht 

We did not bother to introduce new notation for the modified transforms that appear in this theorem 
since these transforms (modified or not) will make no further appearance. We note that the spectral 
transforms introduced in theorem [5] are equivalent to the introduction of coordinates for the bases (in 
the original spectral transform), 

e ±iir/4 / e s/2 \ i 

e A P (s) = t -> ~^== ( Tie -s/2 J S(t, s) fort, s £ T, A P . 
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We follow the developments that lead up to theorem[T] Suppose that s'j 6 £p for j = 1, . . . , m. Define 

ej(s') = ej(s'i) A • • • A ep"(s' m ), where s = (si, . . . , s' m ). 
For Sj € for j = 1, . . . , n define, 

e + A {s) = e+{ Sl )A---Ae+{s n ). 
We apply theorem Q] to the calculation of the multiparticle matrix elements for the spin operator. 

(e+(s'),ae+(s)) =p J R s , xs , R s , xs \ ^ 



(0p,(t0a) V ^ sxs ' R 

where Pf is the Pfafhan and the skew symmetric matrix R has matrix elements, 

_ (e+K)Ae+(*:),<7(U) 

rR v __ (R x _(#>4M_ n -T 

l-Wa'x-Jij- ^xW w - (0 P ,aO A ) ~ <>'i 

and 

Now theorm[T]does not apply because a unitary intertwining map, <r, does not exist as a map from the 
Fock representation on Alt(W^) to the Fock representation on Alt (Wp ). However, all the difficulty 
with the map can be handled by composing it with positive powers of the transfer matrix on both 
sides. The factor (0p,ct0a) remains undefined, however. There are more and less natural ways to 
fix this constant but the discussion makes more sense if one understands how to make a choice for the 
wave function renormalization in the control of the scaling limit. Since we do not want to take up this 
matter here we will set (0p,u0a) = 1 (which is one sensible way in which to scale the spin operator) 
and proceed with the evaluation of (|57|) using the formulas for the matrix elements in theorem [2j We 
begin by defining two (m + n) X (m + n) matrices. Let A denote the (m + n) x (m + n) diagonal 
matrix with non-zero entries, 

A jtj = e~^+ /2 A P (4 + ) for j = 1, . . . , m, 

/ 2mL ch s' 



and 

Next introduce, 
and 



3 



^—mL ch t 



Aj +m ,j+m = 7== e s ' /2 Aa(sj-) for j = 1, . . . 
yj ImL cn Sj 

tj = s'j_ = s[j — in/2, for j = 1, . . . , m, 

t j+m = Sj + = Sj + in/2, for j = 1, . . . , n. 



2G 



Define and (m + n) x (m + ri) matrix, T, by, 

e* 4 — e tj 

Ti j = -7— — 77 fori, j = 1, . . . , m + n. 
Then using the results of theorem [5] we see that, 

Rs'xs' Rs'xs 
Rsxs' Rsxs 

where the transpose A r of the diagonal matrix A is, of course equal to A. Thus we have, 



ATA = ATA T , 



Pf ( p s ' xs ' p s ' xs \ = Pf(ATA T ) = det(A)Pf(T). 

\ -Ksxs' J^sxs J 

See [§J. The determinant of A is just the product of the entries on the diagonal. To finish this 
calculation we observe that the Pfafnan of T has a product form, 

m+n + . + . 

pf ( T ) - n • (58) 

This formula is doubtless well known but can be confirmed in the following manner. Let Z be an 
N x N matrix with entries, 

y _ Z i ~ Z 3 

^i,] ~ — ; — • 
Zi + Zj 

Then Z is a skew-symmetric matrix and, 

N 

mz) = Yl^-^. (59) 

. . Zi -p Zj 

Expand the left hand side by "minors" of the first row (there is a suitable expansion for Pfaffians [8]). 
Use this to determine the residues at the poles (— z^ for k ^ j) of the left hand side in the variable z\. 
We lose nothing by supposing to start that Zi ^ zj fori ^ j; this makes all the poles simple. These 
residues are Pfaffians that can be inductively evaluated using (|59| . Now compare these residues with 
the residues of the right hand side calculated directly. This inductive agreement needs only a check of 
the case N = 2 to be complete. Thus the right and left sides of (|59|) are rational functions with equal 
residues at their simple poles. This implies that the difference of the two sides is a polynomial in Z\. 
It is clear that both sides have finite limits as z\ — > oo and so this polynomial must be a constant. 
Since both sides are at z\ = Z2 they must be the same for all z\. Anytime Zi = Zj for i ^ j, both 
sides of (|59|) are so the equality is true quite generally. 

We are now prepared to state the principal theorem of this paper. 

Theorem 3 Write s' — (s[, . . . , s' m ) with s'j e H' p and s = (si, . . . , s n ) with Sj G T/ A . Recall that for 
^s' =m and #s = n both even, 



and 



e\(s) = ej(si) A • • • A ej(s n ) 
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are both eigenvectors for the transfer matrix. 

The matrix elements of the periodic scaling limit of the Ising spin operator a — ol normalized 
so that (Op,<tOa) = 1 and restricted to a map from A\t CV cn{W A ) to Alt eV en(Wp) in the basis of 
eigenvectors e~p(s') and e A (s) for the transfer matrix is given by, 



'm+n f. i . 

(e+(s'),ae+(s))=A P (s') J[ e -—^-A A (s) 

i<3 



m —imL sh s'j 

A P (s') = e~ im ^ 4 Y[ e- s i/ 2 A P (4 + in/2), 

= i » / 2mL ch s' 



where. 



n p imL sh Sj 

A - ( «>= e,,, " / \n 71sra ^^-' , ^(«»-»/2), 

and 

tj = s'j_ = s'j — in/2, for j = 1, . . . , m, 

tj+m = Sj+ = Sj + in/2, for j = 1,.. . ,n. 

These matrix elements can be used to calculate correlation functions of the spin operators in the 
scaling limit. For the convergence of the resulting sums to be manifest, however, it is helpful if some 
non-zero power of the transfer matrix appears between any two spin operators. We will sketch the 
issues involved for the scaled two point function on the cylinder. 

Return to the lattice formalism on a finite lattice. The two point function on a [21 + 1) x (2m + 1) 
doubly periodic lattice is the ratio of traces, 

_ Tr(<T p V<T V n ) 

{O-p^qO-Qfi) - Tr jy 2 m+1) ' 

where q + n — 2m + 1. Sending m — > oo, fixing p and q one obtains the correlations in the cylindrical 
limit. In this limit the trace in the denominator is asymptotic to A^ m+1 where \a is the largest 
eigenvalue of the transfer matrix. We briefly resurrect the notation from (|3|); don't confuse this with 
the notation A^(s) for the solution to the factorization problem. In the limit n — > oo the product X^ n V n 
tends to the orthogonal projection onto the eigenvector, 0,4, associated with the largest eigenvalue of 
V. Thus in the cylindrical limit the correlation tends to, 

(o-p,qO-o,o) = X A q {0A,(y p V q a Q A )- 

If we introduce a normalized transfer matrix Vp (in the periodic sector) so that Vp — XpVp (and the 
eigenvalue for Vp associated with Op is 1) then the cylindrical correlation becomes, 

(o- p , q o- ,o) = (Xp/\ A ) q (0A,<J P V^ o A ). (60) 

On the finite lattice the spin operator, o~q, is unitary and (Tg = 1. It follows that as long as we work 
in orthonormal bases for Alt CV on(W^4 ) and Alt eV en(M / p~ ) the matrix of o~o thought of as a map from 
Alt GV en(W r p ) to Alt eV en(W^~) is the hermitian conjugate of the matrix of ctq thought of as a map from 
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AltevonCW 7 ^) to Alt eV en(Wp). Replacing cr p by a p = a p /(0 Al a p P ) and cr by ct = cro/(0p, a A ) 
in preparation for scaling, it is straightforward to understand that, apart from space translations, the 
matrices for a p and <to are hermitian conjugates (observe that, modulo the convergence of the scaling 
limit, this shows that it is enough to understand the matrix of the spin operator from Alt eV ea(W^") 
to Alteven(Wp) which is what we've focused on in this paper.) We can now use (f60|) to informally 
understand what happens in the scaling limit to (cr P! q<7o,o)- Make the substitutions q <— qrn^ 1 and 



p <— pm 1 



Note that in the limit we are interested in, 

(x P /x A ) qm ^ ->e-a. 

Thus at least informally, 



pm 1 ,qm 2 



where the sum is over all s — (s±, S2, ■ ■ ■ , s n ) with Sj £ Y,' p and Sj =^ Sj for i =^ j (including s 
Space translation and vertical transfer are given by, 



= J^J exp(imsh(sj)), 



A(s) = Y[ ex P( — mch(sj)). 

The phase factors z(s) p are the only surviving ones in the calculation and one can write, 

2 



\(e+(s),aO A )\ 2 = 



n 



e-^/ 2 A P (s i + m/2) -A- e*- 
11 7*7 



yj/2mL ch 



Results of this sort were announced in [TU] for the scaling functions and derived in [T] for correlations 
on the finite periodic lattice. 
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